Nonlinear response and discrete breather excitation in driven micro-mechanical 

cantilever arrays. 
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We explain the origin of the generation of discrete breathers (DBs) in experiments on damped 
and driven micromechanical cantilever arrays (M.Sato et al. Phys. Rev. Lett. 90, 044102, 2003). 
Using the concept of the nonlinear response manifold (NLRM) we provide a systematic way to find 
the optimal parameter regime in damped and driven lattices where DBs exist. Our results show 
that DBs appear via a new instability of the NLRM different from the anticipated modulational 
instability (MI) known for conservative systems. We present several ways of exciting DBs, and 
compare also to experimental studies of exciting and destroying DBs in antiferromagnetic layered 
systems. 
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The existence and the properties of intrinsic localised 
modes (ILMs) or discrete breathers (DBs) in nonlinear 
lattices have been investigated thoroughly during the 
last years (see p| and references therein). In addition 
an impressive number of various experiments during the 
last years have verified the existence of these modes 
in many systems like micro-mechanical cantilever arrays 
0,0, antiferromagnets 0,0; Josephson junction arrays 
coupled optical waveguides M , atomic vibrations of 
highly nonlinear materials 0, llll and Bose- Einstein 
condensates on optical lattices |l2| . 

The bulk of the central theoretical results has been 
achieved for conservative systems. One good reason for 
that is the complexity of the DB properties. While DBs 
typically persist under the influence of weak dissipation 
(which should also include an energy pumping mech- 
anism), various realisations of dissipation (dc driving, 
ac driving, fluctuations, linear versus nonlinear damp- 
ing etc) modify DB properties in a specific way, turning 
the limiting conservative case into an ideal starting play- 
ground for setting a coherent frame of the understand- 
ing of their properties. Since most of the experimental 
studies face dissipation, each case may call for a specific 
additional theoretical study. 

DB observation in antiferromagnets and can- 

tilever arrays 0, d involved the excitation of the system 
with spatially homogeneous external fields, triggering a 
spatially inhomogeneous system state via some inherent 
instability. For conservative systems the modulational 
instability (MI) of band edge plane waves is known to 
provide such a path. Especially for the case of driven 
cantilevers, the MI approach for conservative systems 
was used to design an experimental system of alternating 
short and long cantilevers 0] • The results we present 
below show that DBs appear via a completely different 
instability. The MI can be used for energy pumping and 
consequent overcoming of dynamical barriers. However 
several other routes with a similar outcome can be ex- 
ploited. We demonstrate that by exciting DBs in arrays 
with identical cantilever length. 

In order to study the driven cantilever system we use 
the model equations in 0, 0] and introduce the dimen- 



sionless time t — > t/to and displacement xi — > xi/ao 
where to and ao correspond to a characteristic time and 
length of the system. The equations of motion describing 
the cantilever system can be then written as a system of 
coupled anharmonic damped and ac driven oscillators 

xi +"/xi +a 2 x t +a±xf —C(xi + i +xi-i —2xi) = A(t) (1) 

The oscillator displacements xi describe the deflection 
angle of the /-th cantilever from its equilibrium position. 
The hard-type anharmonicity tends to increase the oscil- 
lation frequencies with growing amplitudes. This model 
neglects the influence of longer than nearest neighbour 
interaction range, which is not crucial for the under- 
standing of the main qualitative DB properties. The 
dimensionless parameters are related with the ones of 
the experiment in ref.0, 0: 7 = to/r, a i — fato/m, 
04 = k^a^/m, C = ki&jm and A = atg/ao- Using the 
experimental values in [2], |3( and setting a 2 = 0,4 = 1 , we 
find t = 1.34238 • lCT e sec and a = 2.46 • l(T 5 m. The 
friction and coupling parameters become 7 = 1.534- 1CU 4 
and C = 0.07953. The spatially uniform ac driving A(t) 
in m is generated by a corresponding piezoelectric crys- 
tal vibration in the original experiments. 

Neglecting the damping 7, the ac driving A(t) and 
the nonlinear force terms in Q one readily derives 
the only possible solutions, namely plane waves xi ~ 
c i{u q t-qi) wn ^h the linear dispersion relation uj q = 

1 + 4C sin 2 (q/2) relating the plane wave frequency io q 
to its corresponding wave number q. Reinstating the non- 
linear force terms in Q leads to two conclusions i) 
discrete breathers, i.e. time-periodic and spatially lo- 
calised solutions xi{t) — xi(t + T},) , x^^^ — > exist 
for frequencies = > oj q=7r ; ii) the q — it plane 
wave mode turns (modulationally) unstable at ampli- 
tudes ~ 1/N where N is the number of cantilevers, and 
DBs bifurcate from this plane wave mode along this very 
instability route. What can we expect if both damping 
and ac driving are added as well? DBs persist with- 
out much change, thus the ac driving frequency choice 
LJd > u q=T r is well reasoned. However instabilities of con- 
servative systems may turn sensitive to effects of damp- 
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ing. Assuming that the MI of the staggered q = ir mode 
is the track to follow, the choice of a nonstaggered ac 
driving is not appropriate. Note that the experimen- 
tal design of alternating short and long cantilevers splits 
the spectrum Lo q into two bands separated by a gap. In 
addition it formally transforms the q — -k mode of the 
system with identical cantilevers into the q = mode 
of the upper band for the system with alternating can- 
tilever length. Yet the staggered character of this mode 
is of course preserved inside each unit cell which contains 
now two neighbouring cantilevers. Thus the mismatch 
between the staggered MI mode and the nonstaggered 
ac driving remains also for systems with alternating can- 
tilever length. 

With the above parameters 1 < cu q < 1.1481. The 
narrowness of the band as compared to the characteristic 
frequency values is due to the weak coupling constant C. 
Let us start then from the uncoupled limit. In that case 
each oscillator evolves independently. Assuming for the 
moment small amplitudes of all oscillators, the presence 
of a weak driving will (after some transient time) bring 
them all into a unique oscillatory state, thus all oscillators 
will move in phase (with each other). Consequently we 
have to study the stability of a nonstaggered extended 
time-periodic state for weak coupling as well, no matter 
whether the frequency of the driving is located above or 
below the band iv q . 

For a periodic driving of the form A(t) = Aq cos(wdi), 
the equation support periodic solutions. It is easier to 
study first the properties of these solutions at zero fric- 
tion (7 = 0), and then examine the modifications when 
7^0. Newton method 0, ^| was used for the tracking 
of periodic solutions with frequency equal to the external 
driving of the form: xi(t) = xq^i ■ f(u>dt) where xqj is the 
amplitude of the oscillation, and f(ujdt) is a periodic func- 
tion with period 2-7T. The Newton scheme was initiated 
with a very small driving amplitude (Aq ~ 0), and the 
response of the system (i.e. the amplitude of the oscilla- 
tions xqj) was followed as Aq was varying. Thus we re- 
construct the full Nonlinear Response Manifold (NLRM) 

The NLRM close to the origin follows from linearising 
the equation of motion (Q. The exact solution of the 
system in this limit is xi{t) = xq t i cos(tu dt) with xqa — 
—A /(uJ d — 1). This is a homogeneous branch (HB) (all 
the oscillators are in phase). There is a phase difference 
of cj) = 7T between the driving and the response of the 
system. The NLRM is symmetric around the origin due 
to xi(t + ir/uJd) — —Xi(t) and A(t + ir/uJd) — —A(t). 

Increasing \Aq\, the displacement xo.o of HB increase, 
up to the first turning point (TP) (Fig^l. After the 
turning point £0,0 continue to increase, while the driving 
amplitude decreases down to zero at the first crossing 
point (CP1). This crossing point corresponds to a ho- 
mogenous solution with all the oscillators are oscillating 
with the same amplitude xq^q an d zero driving. This state 
exists due to the frequency increasing, hard anharmonic- 
ity term in the equations of motion and because u>d > uj q . 
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FIG. 1: The nonlinear response manifold (NLRM) for C = 
0.0079 and uid = 1.3. Solid line - homogeneous branch, dashed 
line - various breather branches for xo,o- Thick lines corre- 
spond to stable parts of the manifold, thin - to unstable ones. 



The NLRM manifold of the HB continues further from 
CP1 with again nonzero Aq, but the solution is now in 
phase with the driving. 

The Floquet stability analysis [l| of the HB reveals 
that close to the origin the manifold is linearly stable. 
An instability appear before the TP and the HB becomes 
unstable. The NLRM of the homogeneous branch contin- 
ues to be unstable between the TP and the CP1. After 
the CP1, the homogeneous branch turns stable again (see 
FigDJ. 
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FIG. 2: (a,b) The NLRM for different values of the driving 
frequency: (a) ujd = 1.2, (b) uid = 1.4. The parameters are 
the same used in Fig^ ( c ,d) the response £0,0 and the phase 
difference (j> of the manifold in the neighbourhood of the first 
crossing point for different values of the friction. Continuous 
line corresponds to 7 = 0, dotted line to 7 = 0.001 and dashed 
line to 7 = 0.005. 

The instabilities of the HB mark the bifurcation of spa- 
tially inhomogeneous solutions - various breather states 
- off the HB. For a system of TV oscillators, at the first 
instability, N new branches of the NLRM appear. Each 
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of these branches corresponds to a single breather cen- 
tred at a different site on the lattice. In Fig|T]we show 
two of these branches with dashed lines. Starting from 
the bifurcation, the displacement of the central oscillator 
increases, while \Aq\ decreases. For Aq — the manifold 
passes through the second crossing point (CP2). This 
branch of the manifold corresponds to a breather, and is 
stable for Aq > and unstable for Aq < 0. After that, 
the manifold turns, and passes through further crossing 
points. Thus for fixed lo& we have multistability for small 
enough A Q with a large amplitude and small amplitude 
HB coexisting together with breather states. For large 
enough Aq only the large amplitude HB survives. 

The NLRM depends on the driving frequency cud- In 
Fig|21we show a part of the NLRM for two different values 
of the driving frequency. In Fig|2{a) ujd = 1.2 and in 
FigEJb) ujd = 1-4. With increasing frequency the slope 
of the manifold decreases, the crossing points appear for 
larger values of result the bifurcations and the 

turning points appear for larger values of \Aq\. 

The properties of the manifold are slightly modified 
due to the presence of friction (7^0). The NLRM be- 




FIG. 3: Energy as a function of time. The driving frequency 
is initially ramped from ujd — 1.14 to u>d — 1.4, and then it 
remains fixed (marked by the vertical dashed line). Energy 
pumping is observed during the ramping. At the end of the 
ramping process, breathers are created, and one of them locks 
to the driving and survives. T\ = 2Qn/u>d- 

tween the origin and the crossing point (CP1) for 7 = 
is in anti-phase with the driving, while after the cross- 
ing point the NLRM is in phase with the driving. For 
7 = therefore there is a phase difference of 4> — t in 
the response of the system around a CP. When a small 
friction is introduced in the system, it creates an extra 
small phase difference between the driver and the re- 
sponse. This phase difference between the driving and 
the response increases in the neighbourhood of the CP. 
The result is that the solution at the CP disappears due 
to friction, but there is a smooth transition between the 
two different branches of the manifold. The phase differ- 
ence <j) varies smoothly from n — e_ to + e + where e± 
correspond to the small phase differences created by the 
friction on the large and low amplitude branch respec- 
tively, far from the CP. For weak damping this transition 



occurs very close to the CP. The smooth transition be- 
tween the two branches of the manifold in the neighbour- 
hood of (CP1) as a function of \Aq\ is shown in Fig^c) 
for different values of the friction. The smooth transition 
in the phase difference cf> is shown in Fig|2Jd). Far from 
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FIG. 4: (a) Average energy per oscillator as a function of time 
(solid line, scaled by a factor of 10 for better observation) and 
the energy of the locked oscillator (strongly fluctuating dotted 
line). The parameters and the frequency ramping are like in 
FigCD (b) The frequency ramping scheme. After the vertical 
dashed line the driving frequency remains constant. 

the CP, the properties of the manifold are only slightly 
modified by the nonzero friction. 




FIG. 5: (a) Energy density evolution for a sequence of dif- 
ferent driving (see text for details); (b) Time dependence of 
the DB energy observed in (a); (c) Time dependence of the 
DB frequency observed in (a). The two dashed lines mark the 
frequencies wi and U2- The dashed-dotted line corresponds 
to band edge frequency lu w . 

As follows from Fig|21 the stable breather branch is 
not connected to the stable HB of the NLRM. Thus we 
can exclude an easy way of exciting the system in the HB 
and tuning some parameter (amplitude or frequency of 
the driving) such that we continuously join the breather 
branch. This calls for stochastic excitations of the system 
to enforce a nonzero probability to end up in a breather 
state when initially starting from a HB. Similar to the ex- 
periments 2, 3] we performed a ramping of the frequency 
u>d from 1.14 (inside the band) to 1.4 (above the band) for 
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fixed amplitude A = 0.008 (cf. Fig§Jb)). According to 
Fig0the final state corresponds to the multistability do- 
main of the NLRM. If the system was initially at rest, the 
final state is observed to be the low amplitude HB. How- 
ever, adding some nonzero noise, uniformly distributed 
in the interval (—0.05,0.05), in the initial conditions of 
the cantilever deflection angles, the ramping process ex- 
cites a unstable plane wave, which provides a finite time 
energy pumping due to MI (FigJSJ) . The MI acts as a 
noise amplifier, which allows for the generation of hot 
spots or breather precursors. After the ramping finishes, 
most of these hot spots decay back into the low ampli- 
tude HB, while some of them may lock to the driver and 
turn into a stable breather. The average energy per os- 
cillator is shown in Fig^Ja) as a function of time. During 
the ramping, there is a frequency window, where energy 
is pumped into the system, and then the system relaxes 
due to dissipation. The dotted line shows the energy of 
the locked breather site. 

In order to verify that the energy pumping mediated 
by the MI in the experiments is only needed to amplify 
the noise, we perform simulations in the absence of a 
frequency ramping, but with the initial fluctuations be- 
ing larger (Gaussian distribution of initial displacements 
with zero mean and variance D = 0.5). Keeping the 
frequency and the driving amplitude fixed (ujd = 1.4 
and Aq = 0.008) we observe the creation and locking 
of breathers similar to Fig|5] Similar simulations with 
alternating long and short cantilevers, reveal the same 
behaviour. 

The NLRM study shows, that breathers can be ob- 
tained in a driven and damped system by carefully choos- 
ing the frequency (outside the phonon band) and ampli- 
tude of the driving force (not being too large) such that 
the NLRM is activated inside the multistability domain. 
We have shown that localised breathers emerge from the 
homogeneous solution via an instability of the NLRM 
completely different from the expected MI picture. The 
fact that the stable breather branch is disconnected from 



the stable HB, implies that fluctuations have to be used 
in order to perform a crossover from one to the other. 
There exist various pathways of generating breathers, ei- 
ther by frequency ramping and exploiting the MI, or by 
initially strongly exciting the system. These pathways 
have to be designed in such a way that the system will 
tend to the small amplitude HB of the NLRM, with some 
fluctuations growing into hot spots which finally trans- 
form into breather states. 

The NLRM can be also used to study the breather 
excitation in antiferromagnetic systems 0, where 
external driving was used as well. The observed breather 
destruction in these experiments is similar to the process 
presented in FiglHl In analogy with the experiment a 
strong driving at frequency lo\ = 1.8 with A\ = 0.8 
and duration of t = 20Ti is used, together with initial 
displacements being uniformly distributed in the interval 
(-0.25,0.25). As a result a DB is formed (Fi^a)). 
A subsequent long driving at a frequency u>% = 1.4 
(being located closer to the spectrum ui q ) is exciting the 
system with a small amplitude A2 = 0.008. Because 
of the frequency mismatch the breather with lo\ will 
start to decay, but tend to slow down the decay when 
its frequency is close to 0J2 (FhJStc)). Similar to the 
experimental case we add a third weak driving signal 
with 0J3 = 2 and amplitude ^3 = 0.0004. That hinders 
the locking of the DB to oj2- Right after the DB 
frequency passes u>2, its relaxation speeds up and the 
excitation is destroyed very quickly (FigJSJTb), cf. also 
Fig. 3 in £|. The observed energy release is fixed by 
the energy value of a DB locked to the 0J2 driving and 
explains the experimental observation of equal height 
steps in the relaxation of DBs. At the same time we 
interprete the observed smooth decrease of the emission 
signal in Fig. 3 in as a process of slow relaxation of 
DBs with frequency uo\ towards DBs with frequency uj2- 
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